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The gapless color-flavor locked (gCFL) phase is the second-densest phase of matter in the QCD 
phase diagram, making it a plausible constituent of the core of neutron stars. We show that even a 
relatively small region of gCFL matter in a star will dominate both the heat capacity CV and the 
heat loss by neutrino emission L„. The gCFL phase is characterized by an unusual quasiparticle 
dispersion relation that makes both its specific heat cv and its neutrino emissivity e v parametrically 
larger than in any other phase of nuclear or quark matter. During the epoch in which the cooling 
of the star is dominated by direct Urea neutrino emission, the presence of a gCFL region does 
not strongly alter the cooling history because the enhancements of Cv and L„ cancel against each 
other. At late times, however, the cooling is dominated by photon emission from the surface, so 
L v is irrelevant, and the anomalously large heat capacity of the gCFL region keeps the star warm. 
The temperature drops with time as T ~ t~ 1A rather than the canonical T ~ t~ 5 . This provides a 
unique and potentially observable signature of gCFL quark matter. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

It has often been suggested that the core of a neu- 
tron star may contain quark matter in one of the color- 
superconducting phases 0] . The densest predicted phase 
on the QCD phase diagram is the color-flavor locked 
(CFL) phase, which is a color superconductor but an elec- 
tromagnetic insulator 0, Q • The second-densest phase 
is the gapless CFL (gCFL) phase, which is a conductor 
with a nonzero density of electrons [3 It also has 
gapless quark quasiparticles, one of which has an almost- 
quadratic dispersion relation, arising without fine-tuning 
because it is enforced by the requirement that the matter 
be electrically neutral 0, H| ■ We show in this paper that 
this characteristic feature of the gCFL phase means that 
if quark matter in this phase is present in a neutron star, 
it dominates the heat capacity and neutrino luminosity, 
and therefore controls the cooling of the star. At late 
times this produces a unique signature, as the large heat 
capacity keeps the star anomalously warm. A neutron 
star that is tens of millions of years old will be an order 
of magnitude or more warmer if it contains a region of 
gCFL quark matter than if it does not. 

In any speculation about the phases of matter that 
occur inside a neutron star, the main challenge is to pro- 
vide observable signatures of the presence of these phases. 
Since we are proposing such a signature, albeit one that 
presents significant observational challenges, we first set 
the stage with a quick survey of previous proposals. 

• Mass-radius relation. If we could measure the mass and 
radius of several neutron stars to a reasonable accuracy, 
mapping out the mass-radius relationship, we would have 
a strong constraint on the equation of state of dense mat- 
ter. However, although such measurements would dra- 
matically reduce the current uncertainties in our knowl- 



edge of the equation of state of the nuclear matter "man- 
tle" of neutron stars, and could yield evidence of some 
sort of exotic phase in the core, they would not provide 
specific evidence of the presence of quark matter 0, @ • 

• Double pulsar timing. There is a good prospect that 
the long term analysis of the recently discovered binary 
double pulsar [i| may yield a measurement of the moment 
of inertia of a neutron star |ldl | . This would provide infor- 
mation about the density profile that is complementary 
to that obtained from a mass-radius relation, as it would 
constrain the "compactness" of a star. 

• Gravitational waves from collisions. If we could de- 
tect gravity waves from neutron stars spiraling into black 
holes in binary systems, we could perhaps analyze them 
for information about the density profile of the neutron 
star, in particular the presence of an interface separating 
a denser quark core from a less-dense nuclear mantle . 

• Spinning out a quark matter core. If conditions are 
"just so" , rapidly spinning oblate neutron stars may not 
have quark matter in their cores even though more slowly 
rotating spherical neutron stars do. This could be de- 
tected either by anomalies in braking indices of stars 
that are "spinning down" [13 or by anomalous popu- 
lation statistics of stars that are being "spun up" by ac- 
cretion 01 ■ Recent observations show no sign of such an 
effect in the histogram of spin-frequencies of stars in the 
act of being spun up 0| , indicating that if quark matter 
is present, spinning the star and making it oblate does 
not get rid of it. If there is a quark matter core, it must 
therefore occupy a reasonable fraction of the star. 

• r-modes. A rapidly spinning neutron star will quickly 
slow down if it is unstable with respect to bulk flows 
known as r-modes, which transfer the star's angular mo- 
mentum into gravitational radiation. This phenomenon 
will only occur if damping is sufficiently small, so it pro- 
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vides a probe of the viscosity of the interior of the star. 
Such arguments have been used to rule out the possibility 
that pulsars are made entirely of CFL quark matter |l5j , 
in which viscous damping is neglig ible lialig, but their 
implications for the possibility of CFL quark matter lo- 
calized within the core of a neutron star have not yet been 
analyzed. Since gCFL quark matter is expected to have 
a large viscosity, its presence is unlikely to be constrained 
by r-modc arguments. 

• Core glitches. If the third-densest phase on the QCD 
phase diagram is not nuclear matter, it must be a form of 
quark matter with less pairing than in the gCFL phase. 
A leading candidate is the crystalline color superconduct- 
ing phase [T^ j . If this form of quark matter occurs within 
the core of a neutron star, because it is both superfluid 
and crystalline it may be a locus for pulsar glitches. This 
proposal has not yet been worked out sufficiently quanti- 
tatively to determine whether such core glitches are ob- 
servable and if so whether they are consistent with (some) 
observed glitches. 

• Direct neutrino detection. Neutrinos have a long mean 
free path even in nuclear matter, so they can potentially 
carry information about the core directly to the outside 
world. Not coincidentally, neutrinos are very hard to de- 
tect, and the only time when a neutron star emits enough 
neutrinos to be detectable on earth is during the first 
few seconds after the supernova explosion. The time-of- 
arrival distribution of supernova neutrinos could teach 
us about possible phase transitions to and in quark mat- 
ter 0, 0| , but analysis of this proposal requires a bet- 
ter understanding of both the supernova itself and of the 
properties of quark matter at MeV temperatures, where 
the phase diagram of QCD is more baroque than at zero 
temperature ja, @- 

• Cooling. A much better prospect is the indirect de- 
tection of neutrino emission, which is the dominant heat 
loss mechanism for the first million years or so, and can 
therefore be inferred from measurements of neutron star 
temperature as a function of age. Moreover, because both 
neutrino emission rates and heat capacity generally rise 
with density, neutron star cooling is likely to be prefer- 
entially sensitive to the properties of matter in the core 
of a neutron star. 

The qualitative distinction among cooling behaviors that 
may be discerned from the measurement of tempera- 
tures of stars that are 10 3-6 years old is between stars 
in which direct Urea processes are allowed (which yields 
a neutrino emissivity e„ ~ T 6 ), leading to rapid cool- 
ing, and stars in which direct Urea processes are forbid- 
den 0, IH El HI HI HI Ordinary nuclear matter is un- 
usual in that its direct Urea processes n — * p + e + v and 
p + e — > n + v are kincmatically forbidden, meaning that 
neutrino emission relies upon slower reactions (e„ ~ T 8 ). 
Direct Urea processes are allowed in sufficiently dense 
nuclear matter, nuclear matter with nonzero hyperon 
density [2^|, pion condensation (2(| or kaon condensa- 
tion [23, and in all proposed phases of quark matter 
except CFL j2^|. In the CFL phase, there are no direct 



Urea processes because thermally excited quark quasipar- 
ticles are exponentially rare. There are neutrino emis- 
sion processes involving collective excitations that lead 
to e v ~ T 15 0, but in reality any CFL quark matter 
within a star will cool by conduction, not by neutrino 
emission (2^. Indeed, because all forms of dense matter 
are good heat conductors the cooling of a star tends to be 
dominated by whichever phase has the highest neutrino 
emissivity. Hence, the discovery of fast cooling would 
only tell us that some part of the star consists of one of 
the many phases that allow direct Urea. Discovery of 
stars that cool slowly would be an indication that they 
contain only medium-density nuclear matter and perhaps 
CFL quark matter. 



To date, none of the schemes listed above has pro- 
vided an unambiguous signature of the presence of quark 
matter, although all are the subject of ongoing observa- 
tional effort, which in turn drives improvements on the 
theoretical side. In this paper, we argue that recent the- 
oretical advances in our understanding of the properties 
of quark matter offer the prospect of an unambiguous 
detection, if it is possible to measure the temperatures 
of neutron stars that are old enough that their cooling is 
no longer dominated by neutrino emission. Admittedly, 
this presents an observational challenge. However, it is a 
challenge that has not been closely studied prior to our 
work, since all forms of dense matter except gCFL quark 
matter result in neutron stars that cool comparably (and 
very) rapidly in their old age. We show that quark matter 
in the gCFL phase keeps aged neutron stars (those sig- 
nificantly older than a million years) much warmer than 
is predicted by any other assumed dense matter physics. 

Younger neutron stars containing gCFL quark matter 
have a faster-than-standard direct Urea neutrino emis- 
sivity e v ~ T 5 ' 5 , but this does not lead to faster-than- 
standard-direct-Urca cooling because of the correspond- 
ingly enhanced gCFL specific heat. 

In Section ^ we introduce the relevant properties of 
the gapless CFL phase of quark matter, and in Sections 
IIIII and llVl we present the calculations of its specific heat 
cy and neutrino emissivity e v , respectively. These are 
our central calculational results. We do not provide a 
state-of-the-art calculation of the cooling of a neutron 
star containing gCFL quark matter. Instead, in Section 
IV! we provide an introduction to the physics of neutron 
star cooling that suffices to illustrate the qualitative con- 
sequences of the quantitative results for the gCFL cy and 
e u . The astrophysically inclined reader interested in our 
results and their implications but not in their derivation 
can find cy in Eq. (JSJ) and e v in (|31|) and, more con- 
veniently, Fig. |21 and should read the text around these 
results and then turn to Section Ivl 
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II. INTRODUCTION TO THE GAPLESS CFL 
PHASE OF QUARK MATTER 

At any densities that are high enough that nucleons are 
crushed into quark matter, the quark matter that results 
at sufficiently low temperatures is expected to be in one 
of a family of color superconducting phaseSjWith Cooper 
pairing of quarks near their Fermi surfaces [I). The QCD 
quark-quark interaction is strong and is attractive be- 
tween quarks that are antisymmetric in color. If there is 
quark matter in the cores of neutron stars, we therefore 
expect it to be color superconducting. The phenomenon 
persists to asymptotically high densities, where the inter- 
action becomes weak and ah initio calculations of proper- 
ties of color superconducting matter become rigorous [l| . 
The QCD phase diagram exhibits a rich structure of color 
superconducting phases as a function of temperature and 
density QHEl, 

but in this paper we can simplify it by 
working at zero temperature. This is reasonable because 
we will be discussing neutron stars with temperatures in 
the keV range, which is orders of magnitude colder than 
the various critical temperatures at which phase transi- 
tions between different quark matter phases occur. 



A. The CFL phase under stress 

At asymptotically high densities, where the up, down 
and strange quarks can be treated on an equal footing 
and the disruptive effects of the strange quark mass can 
be neglected, quark matter is in the color-flavor locked 
(CFL) phase, in which quarks of all three colors and all 
three flavors form Cooper pairs Q. The CFL phase is a 
color superconductor but is an electromagnetic insulator, 
with zero electron density. In real-world quark matter, 
as may exist in the cores of compact stars, the density 
is not asymptotically high. The quark chemical poten- 
tial \x is of order 500 MeV at most, making it impor- 
tant to include the effects of the strange quark mass M s , 
which is expected to be density dependent, lying some- 
where between the current mass ~ 100 MeV and the 
vacuum constituent quark mass ~ 500 MeV. To describe 
macroscopic regions of quark matter, we must also im- 
pose electromagnetic and color neutrality [3(1 l3ll l3^ | and 
allow for equilibration under the weak interactions. The 
CFL pairing pattern is antisymmetric in flavor, color, and 
spin, so it involves pairing between different flavors. For 
this reason, the effect of a relatively large M s , combined 
with weak equilibration and the neutrality constraints, 
is to put a stress on the CFL pairing pattern: these 
effects would all act to pull apart the Fermi momenta 
of the different flavors by an amount of order j /j, in 
the absence of CFL pairing. (This can be seen by an 
analysis of neutral unpaired quark matter in which the 
Fermi momenta of the d, u and s quarks are split by 
~ Me/4(x.) In the CFL phase, Fermi momenta do not 
separate [33j but the consequence of the stress is that 
the excitation energies of those fcrmionic quasiparticles 



whose excitation would serve to ease the stress by break- 
ing pairs and separating Fermi surfaces is reduced, again 
by of order M^/fi ||. When the density becomes low 
enough, some of the quasiparticles become gapless, the 
CFL pairing pattern is disrupted, and we enter the gap- 
less CFL (gCFL) phase 0,0], the second-densest phase 
on the QCD phase diagram. Since the strength of the 
CFL pairing is measured by the gap parameter Aqfl, 
and the stress on it is of order M s 2 //j, the CFL pattern 
"breaks" and the CFL— >gCFL transition occurs when 
the density is low enough that M s 2 //i ~ Acfl- Making 
this argument quantitative results in the prediction of a 
T = second-order insulator-metal transition separating 
the CFL and gCFL phases at M^/fi ~ 2A C fl 0, H- 
An analogous zero temperature metal insulator transi- 
tion has been analyzed in Ref. l34l. (If the CFL phase is 
augmented by a -RT°-condensate [35|, , the CFL^gCFL 
transition is delaved to a value of Ml / fj, that is higher by 
a factor of 4/3 or less H^-) 

B. The nature of gCFL pairing 

In the gCFL phase the pairing is still antisymmetric 
in flavor as well as color and spin, so as in CFL there 
are diquark condensates (gap parameters) Ai ~ (ds), 
A2 ~ (us), A3 ~ (ud). But unlike the CFL phase, where 
the gap parameters are very similar, Ai = A2 ~ A3, in 
gCFL the pairing involving strange quarks is suppressed: 
very strongly for (ds) , and quite strongly for (us) , so that 
Ai < A 2 < A 3 , as shown in Figures in Refs. 0,0 El- 
The result is that while quarks of all three colors and 
all three flavors still form Cooper pairs, there are regions 
of momentum space in which there is no (ds) pairing, 
and other (very narrow) regions in which there is no (us) 
pairing, and these regions are bounded by momenta at 
which the relevant fermionic quasiparticles are gapless. 

The gCFL phase is an electromagnetic conductor: un- 
like the CFL phase, it contains electrons 0]- The 
electron chemical potential fj, e increases as M^/fi is in- 
creased, rising from zero at the CFL^gCFL transition 
to values which arc comparable to or even larger than 
its typical values in unpaired quark matter, which has 

/i e ~ M S 2 /4 M . 

C. The gCFL domain 

To discuss the range of densities over which gCFL is 
expected to occur, we shall parametrize the strength of 
the attractive interaction between quarks by Ao, which 
we define as the value of the CFL gap parameter at 
M s = in quark matter with fj, = 500 MeV. (We 
shall quote all numerical results at /x = 500 MeV, corre- 
sponding to baryon densities between 8.8 no and 9.1 no 
depending on the value of Ao that we choose, where 
no = 0.17 fm~ 3 is the baryon density in nuclear matter.) 
Because asymptotic-density calculations are not quanti- 



4 




520 540 



p [MeV] 



FIG. 1: Dispersion relations for quasiquarks with gs-bd pair- 
ing (Ai = 3.6 MeV) and 6u-rs pairing (A 2 = 18.5 MeV), in 
the model calculation of Refs. |j. [j.|(J done at /i = 500 MeV, 
M^/u = 100 MeV, with A = 25 MeV. We find gapless 
gs-bd modes at p\ d = 461 MeV and p\ d = 512 MeV. One 
bu-rs mode is gapless with an almost exactly quadratic dis- 
persion relation. Actually it is gapless at two momenta pi u 
and p|", but these are too close together to be resolved until 
the temperature drops below the eV scale, meaning we can 
treat them as a single zero at pi^ = 469 MeV. The five quark 
quasiparticles not plotted are all fully gapped in the CFL and 
gCFL phases. 



gapless, as shown in Fig.^ Wc label the three quark col- 
ors as r,g,b, and make the by now conventional choice for 
which colors pair with which flavors in the CFL phase. In 
this notation one of the gapless branches describes quasi- 
particle excitations that are superpositions of bd particles 
and gs holes. These excitations are gapless at two mo- 
menta p\ d and p^ d shown in Fig. ^ and given by 



Hbd) 



A? 



(1) 



where [i gs and [ibd are determined by the (nontrivial) re- 
quirements of color and electric neutrality. They are de- 
fined in Ref. Q, and their values as a function of M^Ja at 
various Ao can be determined from plots in Refs. [j, |(| . 
The gapless excitations at p\ d are predominantly gs, with 
bd contributing only a small component in the superpo- 
sition. Those at p| arc predominantly bd. In Section llVl 
we shall focus on this dispersion relation in the vicinity 
of p\ d , where it takes the form 



tbd(p) = v bd \p-p bd \ , 



where 



Vbd 



| (pgs - Vbd)] 



(2) 



(3) 



tatively valid at this fi, Aq is not known precisely, with 
estimates ranging from 10 to 100 MeV [|. 

The gCFL phase extends over a range of M^/fi from 
the continuous CFL— >gCFL transition at M| / fi — 2 Ai ~ 
2A up to a first order phase transition where gCFL gives 
way to some phase with even less pairing. In model 
calculations, this transition occurs at M^/fi ~ 5Ao, 
although this is only quantitatively determined within 
particular models |j To give a sense of the 
scales involved, for A = 25 MeV and M s = 250 MeV, 
the gCFL window 2Ao < M^/fi < 5Ao corresponds to 
320 MeV < n < 800 MeV. ' At the lower end of this 
range in [i (upper end in M^/fj,) hadronic matter would 
be more favorable than any form of quark matter. And, 
the upper end of this range in [i (lower end in M 2 //x) cor- 
responds to densities much higher than those achievable 
in neutron stars. Hence, with these choices of parame- 
ters all the quark matter within neutron stars would be 
in the gCFL phase. For larger Ao or smaller M s , the 
gCFL window shifts to lower fi, and neutron stars with 
a CFL core surrounded by a gCFL layer become possi- 
ble. In reality, both Ao and M s are ^-dependent, making 
these estimates illustrative only. 



D. Gapless quasiparticles in the gCFL phase 

In the CFL phase, all nine fermionic quasiparticles are 
gapped. In the gCFL phase, two dispersion relations are 



is the Fermi velocity of the gapless quasiparticles. The 
bd states with momenta between p\ d and p h 2 d arc filled, 
whereas the gs states in this momentum range are empty. 
This means that there is no gs-bd pairing in the ground 
state wave function in this region of momentum space, 
although (as the dispersion relations show) there is still 
pairing among the excitations. The width p| ~~ Pi d °f 
this "blocking region" wherein pairing is "breached" is 
zero at the CFL^gCFL phase transition, and grows 
steadily with increasing M s 2 //x throughout the gCFL 
phase. These dispersion relations behave similarly to 
those describing the gapless modes in the two-flavor gap- 
less 2SC phase H3, and in a metastable three-flavor 
phase discovered in an early analysis in which the con- 
straints imposed by neutrality were not considered |39j . 

The physics of the gapless bu-rs dispersion relation is 
interestingly different. As above these excitations are 
gapless at two momenta p\ u and p h ^ given by 



\iP-bu + Mrs) ± \f [\ {fJ-bu ~ firs)]' 



A 2 



(4) 



but in this instance \\{^bu ~ AVs)] 2 ~ A?, is very small, 
making the dispersion relation in Fig. ^ look quadratic 
with a single zero at p\ u 2 = \{^bu + Mrs)- For the pa- 
rameters of Fig. HI p b 2 u -p^ = 0.026 MeV and the height 
of the dispersion relation half way between these very 
nearby gapless points is only about 5 eV. Since this is 
much smaller than the temperatures that will be of in- 
terest to us, we can safely treat the dispersion relation 
as quadratic, with a dispersion relation in the vicinity of 
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its gapless point given approximately by 



£bu(p) = 



2A, 



(5) 



with the velocity vi> Uj defined analogously to % of 
vanishing at the gapless point. The gap parameter A2 
and the chemical potentials that determine p\ u 2 are plot- 
ted as functions of M 2 //U at several values of A in 
Refs. H@]. 

This near-quadratic dispersion relation is not a result 
of fine tuning. It occurs at all /1 in the gCFL phase, 
and arises from the fact that bulk matter must be elec- 
trically and color neutral 0, 3 . In both the CFL and 
gCFL phase, there is an unbroken gauge symmetry, de- 
noted U{1)q, generated by a linear combination of the 
generators of electromagnetic and color symmetry 0- 
Among the neutrality constraints, it is the imposition 
of Q- neutrality that has the implication of interest. The 
quarks in the gCFL phase are almost Q-neutral by them- 
selves, but not quite: their small excess positive Q charge 
is cancelled by a small admixture of electrons, which have 
Q = — 1. The excess of unpaired M-quarks, occurring in 
a broad band of momenta p bd < p < p 2 d , does not con- 
tribute to the Q imbalance because these quarks have 
Q = 0. It is the unpaired 6u-quarks with p\ u < p < p 2 u 
that matter, because they have Q = +1. They contribute 
a positive Q-charge density of order fj, 2 (p 2 u — p\ u ), bal- 
anced by the electron number density, of order . Since 
[i e *C \i throughout the gCFL phase, p 2 u — p\ u is forced 
(by the dynamics of the gauge fields that maintain neu- 
trality) to remain extremely small, parametrically of or- 
der Pg/^t 2 , throughout the gCFL phase. 

As described above, the dispersion relations for the gs- 
bd quasiparticles are linear about their gapless momenta 
p\ d and p 2 d , as in Eq. J5J), at generic values of M^/fj, in 
the gCFL phase. However, this dispersion relation is fine- 
tuned to be quadratic precisely at the CFL— >gCFL tran- 
sition, where p\ d = p hd and the Fermi velocity Vbd = 0. 
For values of M 2 / /j, that are in the gCFL regime but are 
close to the CFL— >gCFL transition, therefore, the simpli- 
fied linear expression in Eq. @ cannot be used. Indeed, 
if we are interested in those excitations with energies of 
order T or less, the linear expression in Eq. J3J is a good 
approximation as long as vm <] -\/2T/Ai. And, close 
enough to the transition that vm "C y2T/Ai the gs-bd 



dispersion relation can be approximated as quadratic, as 
is appropriate for the bu-rs dispersion relation through- 
out the gCFL phase. 

The gapless excitations of the gCFL phase whose dis- 
persion relations we have described determine the specific 
heat and neutrino emissivity of this phase of matter. In 
the next two sections, we calculate these quantities in 
turn. 

III. SPECIFIC HEAT OF GAPLESS CFL 
QUARK MATTER 

The specific heat of any phase of matter is essentially 
a count of the number of possible excitations with ex- 
citation energies of order T or smaller. Precisely, the 
contribution of quasiparticle excitations with dispersion 
relation e(p) to the specific heat (heat capacity per unit 
volume) at temperature T is given by 



c y = 2 



d 3 p d 
— <P>^ 



(27T) 

2 f d 3 p 

J"2 



dT 



1 



e<P)' T + 1 
<P) 2 



(2tt) 3 ( e <p)/ T + 1) ( e -<P)/ T + l) 



(0) 



where the prcfactor 2 assumes that the quasiparticle is 
doubly degenerate by virtue of its spin. Clearly, only 
those excitations with e < T are important. For the gap- 
less quasiparticle with quadratic dispersion relation, e is 
near zero for p near p\ u 2 and in this regime the dispersion 
relation can be approximated as in Approximating 
the dispersion relation in this way will give us the specific 
heat in the small T limit, and so although it is straight- 
forward to obtain an "exact" result upon assuming JSJ), 
we simply quote the leading result in the small T limit: 



3(V2-1)C(|) 

47r 3/2 

0.146(p^ 2 ) 2 A^ /2 T 1 / 2 . 



1/2 



T l/2 



(7) 



As expected, this is proportional to the number of ex- 
citations with energy less than T, given a quadratic dis- 
persion relation J5j). For quasiparticles with conventional 
linear dispersion relations, namely e(p) = v\p — pp\ for 
some pp and v, the expression JSJ yields the familiar 
cy = jjppT/v. Hence, the specific heat of the gCFL 
phase is given by 



Cy 



{hef 



o.i46 (ply) 2 Al /2 



(k B T) 



1/2 



[Pi 



k B T + 



p bd C y k B T + ^ik B T 



1 ,.2; 



(8) 



where we have restored factors of h, c and Boltzmann's menta are of order /i, and fj-e H- As long as T A 2 , 
constant k B and where the last term comes from the elec- the contribution from the quasiparticle with quadratic 
trons and is negligible because all the quark Fermi mo- dispersion relation dominates. We shall describe reason- 
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able values of A2 and T in subsequent sections; it suffices 
here to say that A 2 /T is of order hundreds or thousands. 

Note that close enough to the CFL— >gCFL transi- 
tion that Vbd < y / 2T/Ai, the gs-bd dispersion relation 
cannnot be treated as linear, and the expression J3J) is 
modified. Indeed, if Vbd <C ^JIT j A\ the gs-bd dispersion 
relation can be treated as quadratic, making their con- 
tribution to the specific heat comparable to that of the 
bu-rs quasiparticles. 

The gCFL phase also has light bosonic excitations, for 
example that associated with superfluidity, but their con- 
tribution to the specific heat is of order T 3 and so can be 
neglected. The specific heat may be enhanced by loga- 
rithimc corrections analogous to those in unpaired quark 
matter |40j |. but we leave their analysis to future work. 



IV. NEUTRINO EMISSIVITY OF GAPLESS 
CFL QUARK MATTER 

The presence of gapless quark quasiparticles in the 
gCFL phase raises the possibility of neutrino emission 
by direct Urea processes. Because the gapless modes are 
superpositions of gs and bd quarks, and bu and rs quarks, 
and because weak interactions cannot change the color of 
a quark, the direct Urea processes that we must consider 



and 



bd — > bu 



bu 



bd 



(9) 



(10) 



Momentum conservation in these reactions requires that 
the momenta of the two quarks and the electron form a 
triangle [2^] . (The argument is that these three fermions 
must all have energies of order T and hence must all be 
close to the momenta at which their dispersion relations 
are gapless. The neutrinos escape from the star, and 
hence have zero Fermi momentum. By energy conser- 
vation, the escaping neutrino therefore has energy, and 
hence momentum, of order T . This is negligible com- 
pared with the momenta of the other fermions.) Mo- 
menta of the two quarks and the electron satisfying the 
triangle constraint can only be found if \p bu — p bd \ < u ei 
where p bu and p bd are the magnitudes of the momenta at 
which the bu and bd quarks are gapless. The only momen- 
tum at which gapless quasiparticles with a bu component 
are found is p\ U 2- Gapless bd quarks occur at two mo- 
menta. For the parameters in Fig. the electron Fermi 
momentum is u e = 25.9 MeV, meaning that the triangle 
constraint can be satisfied if we choose bd quarks with 
momenta near p\ , but cannot be satisfied for those near 
p bd . Indeed, we find that direct Urea processes involving 
the quasiparticles at p bd are forbidden throughout the 
gCFL regime of Mf/11, whereas those involving quasi- 
particles at p bd are allowed throughout all of the gCFL 
regime, with the available phase space vanishing at the 



CFL— >gCFL transition and opening up with increasing 

The calculation of the neutrino emissivity due to direct 
Urea processes in unpaired quark matter was first done 
by Iwamoto in Ref. [28|. and we shall follow his analy- 
sis, leaving the calculation of any logarithmic enhance- 
ment analogous to that in unpaired quark matter |4lj 
to future work. There are two essential differences be- 
tween Iwamoto 's calculation for unpaired quark matter 
and ours for the gCFL phase. First, although the quasi- 
particles near p bd have a conventional linear dispersion 
relation ebdip) given in (J5J), as in unpaired quark matter, 
the quasiparticles near p\ u 2 have a quadratic dispersion 
relation €b u {p) given in (J3J. Analogous to its effect on the 
specific heat, this unusual dispersion relation increases 
the available phase space for the direct Urea reactions by 
a factor of order y/ A2/T relative to that in the standard 
calculation, resulting in a neutrino emissivity e v ~ T 5 5 
rather than ~ T 6 . Second, the quasiparticles near p bd 
and p bu 2 with dispersion relations e^d and tbu are not 
purely bd and bu quarks. They are superpositions of bd 
and gs quarks, and bu and rs quarks, respectively. Only 
the bd and bu components of the quasiparticles partic- 
ipate, meaning that the neutrino emissivity is propor- 
tional to the probabilities that each quasiparticlc is blue. 
These probabilities are the squares of the Bogoliubov co- 
efficients Bb u and Bbd, specified as follows Q- Bb u is 
given by 



1 



BbuipY = o 1 



V 



bu 
Pl,2 



(P 



bu \2 
Pl,2> 



(11) 



where p\ u 2 = \{ub u + Urs) and where A2 is the gap pa- 
rameter for the pairing between bu and rs quarks. To 
simplify the calculation we Taylor expand the coefficient 
around p\ u 2 and get 



Bbu(p) =- 2 + 
The bd Bogoliubov coefficient is given by 



Bbdipf 



1 



1 



p - u 



y/(p-H)*+Al 



112) 



(13) 



where p, = ^(ubd + Hgs) and where Ai is the gap parame- 
ter for the pairing between bd and gs quarks. Because the 
dispersion relation is linear, the quarks that contribute to 
the emissivity lie in a band about p bd whose width is only 
of order T, and we shall replace Bbd{p) by Bbd (Pi d ) ■ This 
coefficient can be quite small. For example, with parame- 
ters as in Fig.Q] meaning in particular / \i = 100 MeV, 
the probability that the gapless quasiparticles at p bd are 

in fact bd is only B bd (p bd ) 2 = 0.00479. 

We now present the calculation of e v , following 
Ref. [2^]. The transition rate for the process © is 



W 



V(27T)^{p bd -p, 



Pbu - p ^ IMI 2 



(14) 
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where the index i runs over the four species that partic- 
ipate in the interaction. V is the normalization volume, 
which will drop out by the end of the calculation, and 
the squared amplitude \M\ 2 is given by 



| Mr = 64G| cos 2 9 c (p bd ■ Pv)( Pbu ■ Pe) 



(15) 



where we have averaged over the spin of the initial down 
quark and summed over the spin of the final up quark. 



Here, Gp is the Fermi constant and 9 C 
angle. 



is the Cabibbo 



The neutrino emissivity is the rate of energy loss per 
unit volume due to neutrino emission. It is obtained by 
multiplying the transition rate by the neutrino energy 
and integrating over the available phase space, weighted 
by the Bogoliubov coefficients. The expression can be 
written as 



J 



2 
V 



.1=1 



d 3 pj 

(2tt) 3 



E v Wn(p bd )(l - n{p bu ))(l - n(p e ))B bu (p bu ) 2 B bd (p[ 



bd\2 



(16) 



r 



Here, the Fermi distribution functions n(p bc i)(l — 
n(jp bu ))(\ — n(p e )) state that in order for the process to 
occur we have to have an occupied down state and un- 
occupied up and electron states. In thermal equilibrium 
(which is maintained by strong and electromagnetic pro- 
cesses occurring on timcscalcs much faster than neutrino 
emission via weak interactions) the distribution functions 
are given by 



i 



exp xi 



where 



where 



Pe - He 

T 



Xbd = ± 



z b d(p, 



Ixl) 



T 



(17) 



(18) 



(19) 



with the ± serving in effect to undo the absolute value 
in J5J), and where 



T 



(20) 



with the ± chosen positive for p bu > p\ u 2 and negative 
for p bu < p\ u 2 . In defining x b d we have used (J5J , meaning 
that this derivation is valid at generic values of M 2 / ji in 
the gCFL regime where v b d > y/2T/A±, but not close 
to the CFL-^gCFL transition, where both the bd and bu 
branches should be treated as quadratic. We shall discuss 
this further below. For later use, we also define 



Pu 

T 



(21) 



It is possible to set the calculation up directly in terms 
of the positive quasiparticle excitation energies, but in- 
troducing the ± as we have done allows us to follow 
Iwamoto's calculation more closely. 

We combine Eqs. JT3J, lfT5jl. ((TrJjl . JTTJl and multiply 
by a factor of 2 in order to include the emissivity due 
to the second process (jlUf) . whose contribution proves to 
be the same as that above. We write the integration 
element d 3 pi = pfdpidVli, where dfli is the infinitesimal 
solid angle. The complete expression for the emissivity 
then takes the form 1281 



C 2 

= TT Acos 2 6» c (l -cose ue )AB . (22) 



16tt 8 



Here, A is an angular integral defined as 



A = (j[ J dn)j S( Pbd - Pbu - p e ) , (23) 



where we have eliminated certain terms that vanish iden- 
tically upon angular integration. A is identical to that 
in Rcf. [23 and upon evaluation yields 



A 



32tt 3 



hd^hu .. 
Pi Pl^r-e 



(24) 



where we have taken |p M | = pf, \p bu \ = p^ 2 , |p e | = fj, e , 
knowing that these are the values at which the B integral 
is dominated. The integral B is defined as 
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poo poo poo poo 

B = Pbddptd Pbudpbu P 2 e dpe I pldp u n(p bd )(l - n(p bu ))(l - n(p e )) 
Jo Jo Jo Jo 

x^<S (xbd — x v - Xbu - x e ) B bu (p bu ) 2 B bd {p bd ) 2 ■ (25) 

Finally, the angle 6 ue in (|22(l is the angle between the momentum of the 6u-quark and that of the electron, when 
the two quark momenta and the electron momentum are arranged in a momentum conserving triangle. A little 
trigonometry shows that 9 ue = 9 de + 9 du where 



{ P \ d ) 2 + £-(jA u 2? 

and 

(p\ d ? + {p\ u 2? 

2p b M 



cos9 de = 7^ZM~. '■ — (26) 



cosd du = TTbdZM ■ ( 2? ) 



Now, all that remains is the evaluation of B. 

The integral B is dominated by p bd near p\ d , by p bu near p\ u 2 , and by p e near fi e so we can pull the factor p 2 d p 2 u p 2 
out of the integrand and replace it by (Pi ) 2 (Pi 2) 2 Me- Next, we change variables of integration from the p's to the 
x's and obtain 



v bd J-oo Jo J-oo y/2T\x bu \ e*™ + 1 e** +le" + l 

x6(x bd + x bu + x e - x v ) (28) 

1/2 

where, as in the calculation of the specific heat, the enhanced density of states dp bu /dx bu = (A 2 T/2 la^l) for the 
quasiparticle with the quadratic dispersion relation is crucial. In (|28[1 we have made the approximation B bu (p bu 2 ) 2 ~ 5, 
since the other term in l)12|) leads to a contribution proportional to T 6 and we are keeping only the leading contribution, 
proportional to T 5 5 . We now rewrite Ij28(l as 

1 7^5.5 D ( n bd\2 roc i-oo roc i-oo 1111 

B = -^={p\ d ) 2 {p\\) 2 A V/X2 * M[Pl ] / dx bd / dx e / dx v xl / dx 



2V _^, v bd J x './ x ./ :: — 7_ 00 -™ ^j^+l^ + le'. + l 

x<5(xb d + £i, u + a; e - !„) , (29) 

use the delta function to perform one of the integrations, and then perform the remaining dimcnsionlcss triple integral 
numerically The result is 

B = 31.18 {pX d ) 2 (A U 2?& V ^ bd[Pl > . (30) 

V b d 

Combining this with the result (|24() for A and substituting into 1221) we obtain the final result for the neutrino 
emissivity of the gCFL phase 

62.36 G\ cos 2 9 C bd bu (k B Tf-^ 2 B bd {p bd f 

£v = 5 TTcT? V 1 ~ COS6» ue )p 1 Pi i2 Me , (31) 

7T° h lv c' v bd 

where we have restored the factors of ft, c and k B - This result is valid as long as v bd > ^/2T/Ai, meaning that the 
gs-bd dispersion relation can be treated as linear with slope v bd . At any nonzero temperature, there is a region just on 
the g'CFL side of the CFL — > gCFL tr ansition where this approximation breaks down. Indeed, in the region so close 
to the transition that v bd <C \/2T/Ai, the gs-bd and bu-rs dispersion relations can both be treated as quadratic, and 
an analysis similar to that we have presented above yields 

Eu = 42 ™ G \ C ° S l 9c (l - cos 9 ue )p\ d p^ 2 MkBTf^A^B bd (p bd ) 2 . (32) 



We shall see below that for temperatures of interest in neutron star physics, this expression is valid only in a 
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M s 2 4i[MeV] 

FIG. 2: Neutrino emissivity e„ of gCFL quark matter from 
igP divided by (r/10 7 K) 5 5 , plotted versus M^/fx. The two 
curves are drawn for two different values of the strength of the 
interaction between quarks, corresponding to CFL gap pa- 
rameters Ao = 25 (solid curve) and Ao = 40 MeV (dashed). 
The location of the CFL — » gCFL transitions for Ao = 25 
and 40 MeV are indicated by dots on the horizontal axis, at 
Ms/fi = 46.8 and 73.0 MeV respectively. To the left of these 
dots, e v is negligible in the CFL phase. As discussed in the 
text we begin the curves a small interval to the right of the 
transition, at the M%/(i where v bd = 0.15. (The A = 25 MeV 
and A = 40 MeV curves begin 2.5 MeV and 3.9 MeV to the 
right of their respective transitions.) The baryon chemical po- 
tential is fi = 500 MeV, corresponding to a density about nine 
times that of ordinary nuclear matter. The effect of changing 
fi while keeping A'f 2 /fi and Ao fixed can be approximated by 
scaling e„ with /i 2 . 

very narrow window of parameter space. It is (|31fl that 
is relevant to neutron star phenomenology. 

The gCFL emissivity (|31|l can be compared to the neu- 
trino emissivity of noninteracting quark matter psj j 

-unpaired 457^ G\ COS 2 9 C 2 6 

- 1680 h"<* M *PF(k B T) 
= (3.6 x 10 14 ergcm _3 s -1 )x 

V 100 MeV / V 500 MeV/ V10 7 K/ 

(33) 

where pf — M is the up quark Fermi momentum and 
where \i e has been replaced by M 2 /4[i, appropriate for 
neutral unpaired quark matter. The gCFL emissivity 
(|31(l is enhanced by a factor of y A2/T relative to that 
of noninteracting quark matter, but the full comparison 
between the two rates is more involved. 

We have obtained our result 1)31(1 in a form which 
makes the dependence of e v on T manifest, but which 
obscures the dependence on M s and [i because A2, 
Bbd{Pi d ), Me, Vbd, V h \ an d Pi"2 all change with M s and 
/j,. The most important fj, dependence is straightforward: 
Pi d ~ Pi U 2 ~ A* and hence e v ~ /1 2 . The remaining depen- 
dence on M s and fi is dominated by the dependence on 



Mg/n, which is nontrivial because A 2 , Bbd(Pi d ), Me and 
Vbd and p^f — p hd all depend nontrivially on M 2 //x. The 
result also depends on Ao, through this same set of quan- 
tities. The reader who wishes to obtain numerical values 
of e v would need numerical values for the gap parame- 
ters and chemical potentials in the gCFL phase given in 
plots in Refs. 0,0, which are used in the specification of 
many of the quantities occurring in i|31|) . Given all these 
implicit dependences, we provide Fig. [21 for the conve- 
nience of the reader who wishes to use our result l|31|) 
for the gCFL neutrino emissivity, for example in order to 
calculate its effects on neutron star cooling. 

The most important dependences of e„ are straight- 
forward: e v ~ T 5 - 5 and e v ~ M 2 - All the remaining 
dependences are best described as dependence on M 2 / fi 
and Ao, and hence are described by Figure [2] which 
shows e„/T 5 ,5 as a function of M 2 /fi for two values of 
Ao- For each Ao, we see nonzero neutrino emissivity 
in the corresponding gCFL regime, with e v negligible at 
lower M 2 //j, in the CFL phase. We do not plot e v at 
values of M 2 //i that are larger than the gCFL regime, 
because it is not known what phase of quark matter 
would be favored there, with what T-depcndence for 
its Ejj. (Crystalline color superconducting quark mat- 
ter is a leading candidate for the third-densest phase 
on the QCD phase diagram, and its neutrino emissivity 
has not been calculated.) Also, at these low densities 
quark matter may well have already been superseded by 
nuclear matter. Note that a reasonable estimate of the 
range of M 2 //J. of interest to describe possible quark mat- 
ter in neutron stars is 50 MeV < M 2 /fi < 250 MeV, 
corresponding roughly to 350 MeV < \i < 500 MeV 
and 150 MeV < M s < 300 MeV. If A is large, say 
A = 100 MeV, the curve on Fig-Oshifts far to the right, 
and any quark matter that occurs is likely CFL, with neg- 
ligible e v . We see from the figure that for Ao = 25 MeV, 
the highest density quark matter that can be reached is 
likely in the gCFL phase. For intermediate values of Ao, 
it is possible to obtain a CFL core surrounded by a gCFL 
layer. 

The shape of the curves in Fig. [21 arises from the com- 
bination of many effects. As M 2 //x increases through the 
gCFL regime, fi e rises monotonically and (1 — cos6 ue ) ini- 
tially rises rapidly as phase space for neutrino emission 
opens up, and then varies little. These effects are over- 
whelmed by the fact that as M 2 / \i increases through the 
gCFL regime, A 2 , 1/vbd, and the Bogoliubov coefficient 
Bbd{p bd ) all decrease monotonically. 

Close to the CFL — > gCFL transition, the most im- 
portant contribution to the steep decrease of £„ with in- 
creasing M 2 //j, seen in Fig.[5|is the factor 1/vbd occurring 
in (13 1 1) . since after all Vbd — at the transition. How- 
ever, one must recall that the expression ED plotted in 
Fig. 2 is only valid for Vbd y^T/Aj". Indeed, for any 
nonzero T there is a region close to the transition where 
Vbd *C y/2T/Ai and the emissivity is given by l|32|) with 
e v - T 5 , not by (ED with e v - T 5 5 as plotted in Fig. [21 
In Fig. [21 we have begun the gCFL curves at the value of 
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Mg/fj, at which vm = 0.15, meaning that the curves can 
be trusted as long as T < Ax/100 ~ A o /100. (Note that 
Ai ~ A near the CFL — > g'CFL transition.) For typ- 
ical neutron star temperatures of order keV, e„ ~ T 5 5 
as given by l|31[l (and the curves of Fig. |3 are therefore 
valid) even closer to the transition than where we stopped 
the curves in the Fig 2. The curves can safely be used 
in neutron star cooling calculations, as we shall do in 
Section El 

Further to the right in Fig. [5J well away from the tran- 
sition, Vbd approaches 1 and the factor 1/vbd ceases to 
control the shape of the curves. In this regime, the most 
important contribution to the decline in e v is the rapidly 
falling Bogoliubov coefficient: as M^/fj, increases p\ d and 
•p hd in Fig. ^ separate and the M-component of the gap- 
less quasiparticle at p\ d , namely Bbd(Pi d ), drops faster 
than ix e rises. 

In SectionEJwe shall sketch the implications of our re- 
sults for the specific heat and neutrino emissivity of g'CFL 
quark matter for neutron star cooling. The most impor- 
tant dependence of e v in this context is its T-dependcncc. 
In all plots in Section we show two curves, both with 
M 2 S I [i = 100 MeV, one with A = 25 MeV and one with 
A = 40 MeV. We choose these values because we see 
from Fig. [21 that they correspond to a small, but reason- 
able, and a larger, but still reasonable, value of e^/T 5 ' 5 . 
Were we to choose values of parameters that happened 
to land very close to the CFL — > gCFL transition, all 
the conclusions that we draw in the next section would 
become stronger. 



V. IMPLICATIONS FOR THE COOLING OF 
NEUTRON STARS 

The central results of this paper are the specific heat 
and neutrino emissivity, calculated in Sections IIIII and 
IIVI We shall not attempt a state-of-the-art neutron star 
cooling calculation here, preferring instead to provide a 
calculation that is better thought of as illustrative, not 
quantitative. The effect that we wish to highlight is both 
large and qualitative, arising from the T-dependence of 
Cy and e v , and we expect that even our crude treatment 
will persuade the reader of its significance. 

We analyze the cooling of a "toy star" consisting of a 
volume of "nuclear matter" at a constant density 1.5no 
and a volume of denser quark matter with a constant 
density specified by /x = 500 MeV. As nuclear matter 
we take an electrically neutral gas of nonintcracting neu- 
trons, protons and electrons in weak equilibrium. We 
investigate three different possibilities for the quark mat- 
ter, all electrically and color neutral and in weak equilib- 
rium, and all with /j, = 500 MeV and M s 2 //i = 100 MeV. 
We consider two possibilities for quark matter in the 
gCFL phase, with A = 40 MeV and A = 25 MeV. 
And, we consider noninteracting quark matter. These 
three options have densities of 9.1, 8.9 and 8.8 times nor- 
mal nuclear matter density no, respectively. By treat- 



ing the quark matter core as having a constant density, 
our calculation neglects the possibility of a thin spheri- 
cal gCFL-CFL interface region, in which there would be 
an enhancement in both the specific heat (by a factor of 
two, see end of Sect. lIIIf) and the neutrino emissivity (see 
Eq. relative to the gCFL expressions |jHJ and l|31|) 

that we shall use. Such a shell would be very thin be- 
cause these enhancements occur only within a very nar- 
row window in Mf/fi, but a quantitative investigation of 
how small its effects are is not possible in our "toy star 
calculation". We choose the quark matter and nuclear 
matter volumes V qm and V nm such that the total mass of 
the star is 1.4 solar masses. If we set the quark matter 
volume to zero, this corresponds to choosing a nuclear 
matter "star" that is a sphere with radius R = 12 km. If 
we include a dense quark matter core with radius i? C ore 
while keeping the total mass fixed, the star shrinks as 
we increase R CO rc- With R corc = 5 km the stellar radius 
is R = 10 km. A gCFL core with radius 5 km has the 
same volume as a gCFL layer extending from r = 4.5 km 
to r = 6 km. Since such a layer would surround a CFL 
quark matter core, and since CFL matter plays no role 
in neutron star cooling, the estimates that we quote for 
^?corc = 5 km can equally well be taken as a guide to 
this scenario. Our final toy-model assumption is that 
our "star" is a black body. The work that needs to be 
done to turn our illustrative "toy star calculation" into a 
quantitative calculation of neutron star cooling includes 
the investigation of realistic density profiles, realistic nu- 
clear matter, and realistic atmospheres. We defer this, 
as our calculation suffices to make our qualitative point. 

Our "star" loses heat by neutrino emission from its 
entire volume and by black body emission of photons 
from its surface. The heat loss due to neutrino emission 
is 

L v = V nm eT + V qm e<T . (34) 

The quark matter neutrino emissivity el" 1 is given either 
by (|31fl or (|33|l , depending on whether we are considering 
a gCFL core or an unpaired quark matter core. The nu- 
clear matter emits neutrinos via modified Urea processes 
like n + X -^p + X + e + is, with X either a neutron or a 
proton that serves to carry away some recoil momentum, 
in order that momentum and energy can both be con- 
served in the process. The resulting emissivity is pH l2fij 

E r = (i.2xio^ Ig c m -3 s -.)(-=.) 3,3 ( T ^) ! . 

(36) 

In evaluating the nuclear matter and quark matter emis- 
sivities, we shall assume that the entire interior of the star 
is at a common temperature T. Both nuclear matter and 
quark matter are good conductors of heat, and neutron 
stars older than a few years are well approximated as 
isothermal. 

Because e v ~ T 5 5 in the g'CFL phase, L„ will be dom- 
inated by neutrino emission from the gCFL matter, un- 
less the gCFL volume is very small. We include cooling 
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curves for cores made of "unpaired quark matter" even 
though this is not expected to be present on the phase 
diagram of QCD at neutron star temperatures because 
it serves as a representative example of the large class of 
phases of dense matter in which e„ ~ T 6 and cy ~ T. 
This class includes all quark and nuclear phases that cool 
by direct Urea processes, except for gCFL. 

The surface of real neutron stars is colder than their 
interiors, with the temperature gradients occurring only 
in the outer envelope of the star within of order 100 me- 
ters of the surface. The heat transport within this en- 
velope has been analyzed ^3]; an d the result is well ap- 
proximated by a phenomenological relationship between 
the interior temperature T and the surface temperature 
T surfacc given by [HIH 



T surfa cc = (0.87 x 10 6 K) 



9s 



10 14 cm/s : 



1/4 



0.55 



(36) 

where g s = GnM/R 2 is the surface gravity. This means 
that the rate of heat loss from the surface of the star, 
which for a black body is 



L 7 = 47ri? 2 crT s 4 urfacc 
with a the Stefan-Boltzmann constant, is given by 



(37) 



L 1 = 4irR 2 a(0.87 x 10 6 K) 4 



fjs 



10 14 cm/s 5 



(—)' 
\ 10 8 K/ 



(38) 

We shall use this expression for L 7 , even though we are 
not treating other aspects of the problem realistically, 
because the fact that L 1 ~ T 2 - 2 will play an important 
qualitative role. 

The cooling of our "star" is described by the differen- 
tial equation 

dT L v + L 1 

df ~ V r nrn 4- V r qm 
u< h v nm^y \ v qm^y 



V F 



nm ,y qm , r 



y „nrn , y qm 

(39) 

which equates the heat lost (~ Ldt) to the change in 
the heat energy of the star (~ — VcydT). We have all 
the ingredients needed to evaluate the right hand side of 
this equation in place, with the exception of the specific 
heat of nuclear matter and of unpaired quark matter. 
For a gas of several species of nointeracting fermions, the 
specific heat is given by 



klT 



(40) 



where the sum runs over all the species. In the case 
of noninteracting nuclear matter, the sum runs over i = 
n, p, e and the Fermi momenta for neutral matter in weak 
equilibrium are given by [2l| 



p n F = (340 MeV) [ — 



\n 



1/3 



P P F = p) = (60 MeV) 



2/3 



(41) 
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FIG. 3: Solutions to the cooling equation for 1.4 solar 
mass "toy stars" (described in the text) of four different com- 
positions. The curves show internal temperature as a func- 
tion of time. The dot-dashed curve is for a star with radius 
R = 12 km made entirely of nuclear matter with a density 
1.5no, with no quark matter core. The other three curves 
describe stars with radii R — 10 km that have quark mat- 
ter cores with radii 7? C orc = 5 km. For all three curves, the 
quark matter has fx = 500 MeV and M 2 S / > = 100 MeV, with 
densities ~ 9no. For the dotted curve, the quark matter is 
noninteracting. For the solid (dashed) curve, it is in the gCFL 
phase with A = 25 MeV (A = 40 MeV). 



We are taking noninteracting nuclear matter with density 
n = 1.5no for the mantle of our "stars". In the case 
of neutral unpaired quark matter in weak equilibrium, 
the sum on i runs over the nine quarks and the Fermi 
momenta are independent of color and are given by |30l | 



Pf 



Pf 
Pf 



Ml 

12yU 

Ml 

6/x 
5Ml 
Yip, 



(42) 



up to corrections of order M '£ / fi 3 . We are using matter 
with /i = 500 MeV in the core of our "stars" . 

Fig-EDshows the cooling curves obtained by solving i|3"9"|) 
for the four toy stars we have described, plotted on a log- 
log plot. Each curve has an early time power law during 
the period when cooling by neutrino emission dominates, 
namely the first 10 5 or so years. At early times, L 7 « 
because L 7 ~ R 2 whereas L v ~ R 3 . Because L v drops 
much more rapidly than L 7 as T decreases, at late times 
.L 7 dominates and a new power law is seen. 

It is easy to sec why power law solutions arise. In 
any temperature regime in which the numerator and the 
denominator of the right hand side of l|39|l are each domi- 
nated by one of their component terms, the cooling equa- 
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occur that has been considered previously. But, the ef- 
fect of these enhancements cancel in the cooling curve 
during the epoch when L„ > L 7 . There are now a num- 
ber of indications that some neutron stars with ages 
of order 10 3 to 10 5 years (presumably the heavier ones, 
although this is certainly not demonstrated) are signifi- 
cantly colder than would be expected in the absence of di- 
rect Urea neutrino emission, whereas other (presumably 
less massive) stars have temperatures consistent with the- 
oretical cooling curves calculated upon assuming nuclear 
matter composition. Were this to be confirmed, the dis- 
covery of direct Urea emission with T ~ t -1 / 4 , instead of 
the slower T ~ t -1 / 6 , could indicate the presence of any 
number of dense matter phases, including gCFL quark 
matter but also including nuclear matter leavened with 
either hyperons, kaons or pions. 



FIG. 4: Same as Fig. [31 except that here we plot T B urfacc, 
related to the interior temperatures plotted in Fig. |H] by 
Eq. 



tion takes the form 

f = ™ 

for some p and a. For example, for a star that is made 
entirely of nuclear matter, during the epoch when L v 3> 
L 7 we have p = 7, coming from L y ~ T 8 and c v m ~ T. 
For p > 1, (|43[) has a power-law solution 

T= [a(p-l)t]~7^ . (44) 

There are no arbitrary constants in this solution. We ini- 
tialize the differential equation with some temperature To 
at a time t = I year, chosen because by that time the in- 
terior star can reasonably be treated as isothermal. The 
initial condition To (to) does not appear in the power law 
solution: it only affects how the power law solution is 
reached, if To (to) does not lie on it. Once the power law 
solution is reached, the form of the solution to the dif- 
ferential equation is independent of the initial condition. 
We begin all our plots at t = 1000 years, by which time 
the solution is on the power law (|44|) for any reasonable 
choice of T (t ). 

During the epoch when L v L 1 , a star made entirely 
of nuclear matter has p — 7 and T ~ t -1 / 6 . For the 
stars with unpaired quark matter, or gCFL quark mat- 
ter, p = 5 and T ~ t -1 / 4 during this epoch. This explains 
how similar the three quark matter core curves are dur- 
ing the first 10 5 years, and why all three stars with quark 
matter cores are colder than the nuclear matter star. 
Note that the gCFL quark matter has L v ~ T 5 5 and 
cy ~ T 5 , both enhanced by l/T 5 relative to that of 
unpaired quark matter, and indeed relative to any phase 
of nuclear or quark matter in which direct Urea processes 



At late times, when L 7 ^> L u all stars except those 
containing gCFL quark matter have p = 2.2 — 1 = 1.2 
because L 1 ~ T 2,2 and cy ~ T, and hence cool with 
T ~ t -1 / - 2 = i~ 5 . This explains the rapidly dropping 
temperatures at late times for the stars without gCFL 
quark matter in Fig. [3J If the volume of gCFL matter 
is sufficient (more on this below) it dominates the heat 
capacity of the star, yielding p = 2.2 — 0.5 = 1.7 because 
Cv ~ T 5 , and hence the star cools with T ~ t" 1 / 7 = 
t~ 1A at late times. The gCFL matter keeps the aging 
star warm by virtue of its large heat capacity. Hence, 
the title of our paper. 

We show the surface temperatures of our toy stars in 
Fig. El It is tempting to put data obtained from the 
observation of real stars on this plot, but we resist the 
temptation given that our "stars" are not realistic. The 
qualitative impact of gCFL quark matter is, however, 
clear: stars which are old enough that they cool by pho- 
ton emission stay much warmer if they contain a gCFL 
hot water bottle. In our Fig.0J which should be taken as 
illustrative and not yet as a quantitative prediction, the 
effect is a full order of magnitude for 10 7 -year old stars, 
and gets much larger for older stars, as the cooling curves 
of all stars except those containing gCFL quark matter 
drop rapidly. 

In Fig. we investigate the dependence of the cooling 
curves on the volume of gCFL quark matter present in 
the core of the star. We see that the "hot water bottle 
effect" is present for i? corc = 3 km, but reduced in mag- 
nitude. For i? corc = 1 km, no effect is visible: the effect 
does occur, but only at even lower temperatures than we 
have plotted. (Because its heat capacity is cy ~ T 5 , if 
any gCFL quark matter is present it will eventually dom- 
inate the heat capacity of the entire star, no matter how 
small its volume fraction. For R COTC = 1 km, this occurs 
at temperatures below those we have plotted.) Note that 
what we are referring to as R COTC = 5 km could equally 
well describe a star with a shell of gCFL quark matter 
extending between radii of 4.5 and 6 km. 
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FIG. 5: Cooling curves showing the surface temperature of 
stars with g'CFL cores with Ao = 25 MeV for i? C orc = 5 km 
(solid; same as solid curve in Fig. 0), -Rcoro = 3 km (dashed), 
Rcoia = 1 km (dotted), -Rcoro = km (dot-dashed; same as 
dot-dashed curve in Fig. . 



VI. OUTLOOK 

We hope that our results challenge observers to con- 
strain the temperature of neutron stars that are 10 mil- 
lion years old or older. Prior to our work, all proposed 
cooling curves for these old stars drop so fast into unob- 
scrvability that there has been little motivation to make 
the effort to obtain the best constraints possible on their 
temperatures. Given that we know of isolated neutron 
stars that are younger than a million years old and closer 
than 200 parsecs, it is reasonable to expect that there are 
10 7 -year old isolated neutron stars closer than 100 par- 
secs to earth. At first we were concerned that if such old 
nearby stars had temperatures of order 10 5 K, as Fig. 0] 
suggests will be the case if they contain gCFL hot water 
bottles, they should already have been detected. How- 
ever, initial estimates suggest that they will in fact be 
quite a challenge to find, since the peak of a 10 5 K black 
body spectrum lies in the far ultraviolet, where the in- 
terstellar medium is opaque, and since they will be quit e 
faint in the accessible UV and visible wavelengths |44j . 
Another option, perhaps easier than finding these stars 
without knowing where to look, is to study nearby old 
pulsars, already detected by their nonthermal emission, 
and to constrain their thermal emission hence putting an 
upper bound on their temperature. This has been done 
for PSR 0950+08, whose spin-down age is 10 7 2 years, 
yielding the bound T < 10 5 2 K 45]. This limit is quite 
promising, as it is close to the curves in Fig. 01 describing 



the cooling of our toy star with a gCFL core. And, we 
are confident that we have not thought of the best way 
of looking for aging but still warm neutron stars. We are 
therefore hopeful that the opportunity to make an unam- 
biguous discovery of the presence of quark matter within 
neutron stars or to rule out the presence of gCFL quark 
matter in the entire region of the QCD phase diagram 
sampled by neutron stars will stimulate observers to rise 
to the challenge. 

Much theoretical work remains to be done. Interesting 
microphysical questions about the gCFL phase remain, 
and have been enumerated in Refs. 0, |(| . Phases with 
some features in common with the gCFL phase can be un- 
stable with respect to inhomogeneous mixed phases |46j| , 
and although the gCFL phase is stable with respect to all 
straightforward mixed phase possibilities 0, w\ , an ex- 
haustive investigation has not yet been performed. Per- 
haps the most interesting open ques tions are the pos- 
sibilities of A'°-condensation pjrj l37l| or gluon conden- 
sation [4^ in the gCFL phase. Either could change our 
quantitative results for its cy and e„, but neither is likely 
to change their unusual T-dcpcndcnce: cy ~ T 5 and 
e v ~ T 5 5 . (Neither if°-mesons nor gluons [3(], H2 would 
affect the Q-chargc balance, which is responsible for the 
existence of the gapless quasiparticle with a quadratic 
dispersion relation in Fig. ^ whose consequence in turn 
is the unusual T-dependence of the gCFL cy and e„.) As 
far as theoretical astrophysical work, our results for cy 
and e v must be incorporated into calculations of cooling 
curves for stars with realistic atmospheres and density 
profiles before plots like those in Figs. 0] and can be 
compared quantitatively to data. Nonetheless, our con- 
clusion that gCFL quark matter within a neutron star 
will keep the star warm in its old age relies only on the 
unusual T-dependence of the gCFL specific heat, and is 
therefore expected to be robust. 



Acknowledgments 

We thank David L. Kaplan for orientation on various 
questions related to observations and Kenji Fukushima 
for assistance with quasiparticle dispersion relations, and 
acknowledge helpful conversations with them and with 
Michael Forbes. We thank Igor Shovkovy for finding 
an error in a previous version of this paper. PJ is 
grateful to the Research Science Institute of the Center 
for Excellence in Education for supporting her research. 
This research was supported in part by DOE grants DE- 
FG02-91ER40628, DF-FC02-94ER40818, and DE-FG02- 
93ER40762. 



[1] For reviews, see K. Rajagopal and F. Wilczek, 
|arXiv:hep ph/0011333| M. G. Alford, Ann. Rev. Nucl . 
Part. Sci. 51, 131 (2001) |arXiv:hep-p h/0102047 ; 



G. Nardulli, Riv. N uovo Cim. 25N3, 1 (2002) 
[arXiv:hep-p h/0202037 ; S. Reddy, Acta Phys . 
Polon. B 33, 4101 (2002) |arXiv:nucl-th/0211045| ; 



14 



[2] 
[3] 
[4] 

[s; 
[(■> 
[r 
[§ 

[9 
[10 

[11 

[12 
[13 
[14 

[15 
[1(3 
[17 



T. Schafer, |arXiv:hep -ph/0304281 M. Alford, 
arXiv:nucl-th/0312007 

M. G. Alford, K. Rajagopal a nd F. Wilczek , Nucl. Phys. 

B 537, 443 (1999) |arXiv:hep- ph /9804403] . 

M. Alford, C. Kouvaris and K. Rajagopal, Phys. Rev. 

Lett. 92, 222001 (2004) arXiv:hep-ph/0311286 . 

M. Alford, C. Kouvaris and K. Rajagopal, Phys. Rev. D., 

to appear arXiv:hep-ph/0406137 . 

S. B. Ruster, 1. A. Shovkovy and D. H. Ris chke, Nucl. 

Phys. A 743, 127 (2004) arXiv:hep-ph/0405170 . 

K. Fukushima, C. Kouvaris and K. Rajagopal, 

arXiv:hep-ph/0408322 

J. M. Lattimer and M. Prakash, Science 304, 536, (2004) 



arXiv:astro-ph/0405262 

m7 



M. Paris and S. Reddy, 



Alford, M. Brab y, 
|arXiv:nucl-th /041 i016] 

M. Burgay et al, Nature 426. 531 (2003); A. G. Lyne et 

ai, Science 303, 1 153 (2004). 

M. Kramer et al, |arXiv:astro- ph/0405179 I. A. Morri- 
son, T. W. Baumgarte, S. L. Shapiro and V. R. Pand- 
haripande, arXiv:astro-ph/0411353 J. M. Lattimer and 

B. F. Schutz , |arXiv: astro-ph / 041 1470| 
M. G. Alford, K. Rajagopal, S. Reddy and F. Wilczek , 
Phys. Rev. D 64, 074017 (2001) |arXiv:hep-p h/0105009 . 
N. K. Glendenning, S. Pei and F. Weber, Phys. Rev. Lett. 
79, 1603 (1997) |arXiv:astro-ph/9705235| . 
N. K. Glendenning and F. Weber, Astrophys. J. 559, 
L119 (2001) |arXiv:astro- ph /0003426] . 
D. Chakrabarty et al., Nature 424, 42 (2003); R. Wi- 
jnands et al, Nature 42 4, 44 (2003); D. Chakrabarty, 
arXiv : astro-ph /0408004 

J. Madsen, Phys. R ev. Lett. 85, 10 (2000) 
arXiv:astro-ph/9912418 . 

C. Manuel, A. Dobado and F. J. Llanes-Estrada, 
arXiv:hep-ph/0406058 

M. G. Alford, J. A. Bowers and K. Rajagopal, Phys. Rev. 
D 63, 074016 (2001) |arXiv:hep-p h /0008208 ; J. .A. Bow- 
ers, J. Kundu, K. Rajagopal and E. Shuster, Phys. Ref. 
D 64, 014024 (2001) larX"iv:hep-ph/0101067| ; J. Kundu 
and K. Rajagopal, Ph ys. Rev. D 65, 094022 (2002) 
arXiv:hep-ph/0112206 ; J. A. Bowers and K. Rajagopal, 
Phys. Rev. D 66, 065002 (2002) arXiv:hep-ph/0204079 ; 
R. Casalbuoni and G. Nardulli, Rev. Mod. Phys. 
263, 320 (2004) arXiv:hep-ph/0305069 . R. Casalbuoni, 
M. Ciminale, M. Mannarelli, G. Nardulli, M. Ruggieri 
and R. Gatto, |arX"iv:hep-ph/ 0404090 
[18] G. W. Carter and S. Reddy~Phys. Rev. D 62, 
103002 (2000) arXiv:hep-ph/0005228 ; S. Reddy, 
M. Sadzikowski a nd M. Tachibana, Nu cl. Phys. A 
714, 337 (2003) arXiv:nucl-th/0203011 ; S. Reddy, 
M. Sadzikowski and Tachibana, Phys. Rev. D 68, 
053010 ( 2003) |arXiv:nuc l-th/0306015 ; J. Kundu and 
S. Reddy, |arX"iv:nucl-th/0405055| 
[19] P. Jaikumar, M. P rakash and T. Schafer, P hys. Rev. D 

66, 063003 (2002) |arXiv:astro~p h/0203088] . 
[20] G. Baym and C. J. Pethick, Ann. Rev. Astron. Astro- 
phys. 17, 415 (1979); C. J. Pethick, Rev. Mod. Phys., 64, 
1133 (1992); M. Prakash, in Nuclear and Particle Astro- 
physics, eds. J. G. Hirsch and D. Page, 153 (Cambridge: 
Cambridge University Press, 1998); D. G. Yakovlev, 
A. D. Kaminkcr, O. Y. Gnedin and P. Haensel, 
Phys. Rept. 354, 1 (2001) arXiv:astro-ph/0012122 ; 

D. G. Yakovlev and C. J. Pethick, Ann. Rev. Astron. 
Astrophys. 42, 169 (2004). 



[21 
[22 

[23 
[24 



[25 
[26 

[27 



[28 
[29 
[30 
[31 



S. L. Shapiro and S. A. Teukolsky, Black Holes, White 
Dwarfs and Neutron Stars, (New York: Wiley, 1983). 
D. Page, M. Prakash, J. M. Lattimer and 
A. Steiner, Phys. R ev. Lett. 85, 2048 (2000) 
[arXiv:hep- ph /0005094] . 

D. Page, J. M. L attimer, M. Prakash and A. W. Steiner, 
|arXiv:astro-ph/0403657l 

D. G. Yakovlev, O. Y. Gnedin, M. E. Gusakov, 
A. D. Kaminker, K. P. Levenfish and A. Y. Potekhin, 
arXiv:astro-ph/0409751 

N. K. Glendenning, Astrophys. .1., 293. 470 (1985). 

J. N. Bahcall and R. A. Wolf, Phys. Rev. B 140, 1452 

(1965). 

D. B. Kaplan and A. E. Nelson, Phys. Lett. B 175, 57 
(1986); D. B. Kaplan and A. E. Nelson, Phys. Lett. B 
179, 409E (1986) 

N. Iwamoto, Phys. Rev. Lett. 44, 1637 (1980); 
N. Iwamoto, Annals Phys. 141, 1 (1982) 
I. A. S hovkovy and P. J. Ellis Phys. Rev. C 66, 015802 
(2002) |arXiv:hep-ph/0204132| . 

M. Alford and K. Raja gopal, JHEP 0206, 031 (2002) 
[arXiv:hep-ph/020400l] ; 

Phys. Rev. D 63, 074018 
A. W. Steiner, S. Reddy 
Rev. D 66, 094007 (2002) 
F. Neumann, M. Buballa 
Phys. A 714, 481 (2003) 
A. Gerhold and A. Rebhan, 



K. lida and G. Baym, 
(2001) |hep-ph /0011229 ; 
and M. Prakash, Phys. 

arXiv:hep-ph/0205201 ; 
and M. Oertel, Nuc l. 

arXiv:hep-ph/0210078 ; 



[32] 
[33] 

[34] 

[35] 



[36] 
[37] 

[38] 



[39] 
[40] 



[41] 
[42] 



[43] 



Phys. Rev. D 68, 011502 (2003) arXiv:hep-ph/0305108 ; 
A. Kryjevski Phys. Rev. E» 68, 074008 (2003) 
|arXiv:hep-ph/0305173| ; D. D. Dietrich and 

D. H. Rischke, Prog. Part Nucl. Phys. 53, 305 
(2004) |arXiv:nucl-t h /0312044] . 

A. Gerhold, |arXiv:hep-ph/04Tl 086 

K. Rajagopal and F. Wilczek, Phys. Rev. Lett. 86, 3492 
(2001) |arXiv:hep-p' h/001203^1. 

W. V. Liu, F. Wilczek and P. Zoller, 
|arXiv:cond-mat/0404478l 

P. 1'. Bedaq ue and 'I\ Schafer, Nu cl. Phys. A 697, 
802 (2002) |arXiv:hep-ph/0105150| ; D. B. Kaplan 
and S. Reddy, Phys Rev. D 65, 054042 (2002) 
[arXiv:hep-p h/01072 6'5]; A. Kryjevski, D. B. Kaplan 
and T. Schafer, |arXrv7h ep-ph/0404290 M. Buballa, 
|arXiv:hep- ph/0410 397[ 

M M. Forbes, |arXiv:hep-ph/041100T1 

A. Kryjevski and T. Schaef er, |arXiv:hep-ph/04 07329 
A. Kryjevski and D. Yamada, arXiv:hep-ph/0407350 
I. Shovkovy and M. Huang, Phys. Lett. B 564, 205 
(2003); E. Gubankova, W. V. Li u and F. Wilczek, Phys . 
Rev. Lett. 91, 032001 (2003) arXiv:hep-ph/0304016 ; 
M. Huang and I. Shovkovy, |arXiv:hep -ph/0307273 
M. G. Alford, J. Berges and K. Rajagopal, Phys. Rev. 
Lett. 84, 598 (2000) arXiv:hep-ph/9908235 . 
A. Ipp, A. Gerhold an d A. Rebhan, Phys. Rev. D 
69, 011901 (2004) |arXiv:hep- ph/03090 19|; A Gerhold, 

A. Ipp and A. Rebhan, |arXiv:hep-ph/0406087| 

T. Schafer and K. Schwenzer, arXiv:astro-ph/0410395 

E. H. Gundmundsson, C. J. Pethick and R. I. Epstein, 
Astrophys. J. 259, L19 (1982); E. H. Gundmundsson, 
C. J. Pethick and R. I. Epstein, Astrophys. J. 272, 286 
(1983). 

P. Slane, D. J. Helfand and S. S. Murray, Astrophys. 
J. 571, L45 (2002) [arXiv:astro-ph/020415lT ; D. L. Ka- 
plan et ai, Astrophys. J. Supp. 153, 269 (2004) 



15 



arXiv:astro-ph/0403313 ; J. P. Halpern, E. V. Gotthelf, 
F. Camilo, D. J. Helfand and S. M. Ransom, Astrophys. 
J. 612, 398 (2004) [astro-ph/0404312| . 

[44] D. L. Kaplan, private communication. 

[45] V. E. Z avlin and G. G. Pavlov Astrophys. J. 616, 452 
(2004) arXiv:astro-ph/0405212 . See D. Page, Astro- 
phys. J. 479, L43 (1997) and Ref. H3| for an attempt 
to obtain such a limit using older data. 

[46] P. F. Bedaque, Nuc l. Phys. A 697, 569 (2002) 
|arXiv:hep-ph/9910247| ; P. F. Bedaque, H. Cal- 
das an d G. Rupak, Phys. Rev . Lett. 91, 247002 
(2003) |arXiv:cond-m at/0306 6'94l; H. Caldas, Phys . 
Rev. A 69, 063602 (2004) |arXiv:hep-ph/0312275| ; 
M. M. Forbes, E. Gubankova, W. V. Liu and 



F. Wilczek, |arXiv:hep-ph/0405059| S. Reddy and G. Ru- 
pak, |ar»7inucn^704050Mp 

[47] M. Alford, C. Kouvaris and K. Rajagopal, 
arXiv:hep-ph/0407257 

[48] M. Huang and 1. A. Shovkovy, Phy s. Rev. D 70, 
051501R (2004) arXiv:hep-ph/0407049 ; M. Huang and 
I. A. Shovkovy, Phys . Rev. D 70, 094030 (2004) 
arXiv:hep-ph/0408268 ; R. Casalbuoni, R. Gatto, 
M. Mannarelli, G. Nardulli and M. Ruggieri, Phys. Lett. 
B 605, 362 (2005) |arXiv:he p-ph/0410401 ; I. Gian- 
nakis and H. C. Ren, Phys. Lett. B 611, 137 (2005) 
|arXiv:hep-ph/0412015| ; M. Alford and Q. h. Wang, 
|arXiv:hep-ph/0501078T 



